Symplectic Molecular Dynamics Integration Using Normal Mode Analysis
نویسندگان
چکیده
The split integration symplectic method (SISM) for molecular dynamics (MD) integration using normal mode analysis based on a factorization of the Liouville propagator is presented. This approach is quite distinct from others that use fractional-step methods, owing to the analytical treatment of high-frequency motions. The method involves splitting the total Hamiltonian of the system into a harmonic part and the remaining part. Then the Hamilton equations are solved using a second-order generalized leapfrog integration scheme in which the purely harmonic Hamiltonian (which represents the main contribution of the chemical bonds and angles) is treated analytically, i.e., independent of the step size of integration, by a normal mode analysis that is carried out only once, at the beginning of calculation. The whole integration step combines analytical evolution of the harmonic part of the Hamiltonian with a correction arising from the remaining part. The proposed algorithm requires only one force evaluation per integration step. The algorithm was tested on a simple system of linear chain molecules. Results demonstrate the method makes possible the integration of the MD equations over larger time steps without loss of stability while being economical in computer time. c © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 2–12, 2001
منابع مشابه
Molecular dynamics integration and molecular vibrational theory. I. New symplectic integrators.
New symplectic integrators have been developed by combining molecular dynamics integration with the standard theory of molecular vibrations to solve the Hamiltonian equations of motion. The presented integrators analytically resolve the internal high-frequency molecular vibrations by introducing a translating and rotating internal coordinate system of a molecule and calculating normal modes of ...
متن کاملMolecular Dynamics Integration Meets Standard Theory of Molecular Vibrations
An iterative SISM (split integration symplectic method) for molecular dynamics (MD) integration is described. This work explores an alternative for the internal coordinate system prediction in the SISM introduced by JaneZic et al. (J. Chem. Phys. 2005, 122, 174101). The SISM, which employs a standard theory of molecular vibrations, analytically resolves the internal high-frequency molecular vib...
متن کاملParallel Approaches in Molecular Dynamics Simulations
In this contribution we will present the survey of our past and current endeavor on parallel approaches in molecular modeling algorithm development, for example, molecular dynamics (MD) simulation. In particular, we will describe the new split integration symplectic method for the numerical solution of molecular dynamics equations and methods for the determination of vibrational frequencies and...
متن کاملPreservation of Adiabatic Invariants under Symplectic Discretization
Symplectic methods, like the Verlet method, are a standard tool for the long term integration of Hamiltonian systems as they arise, for example, in molecular dynamics. One of the reasons for the popularity of symplectic methods is the conservation of energy over very long periods of time up to small fluctuations that scale with the order of the method. In this paper, we discuss a qualitative fe...
متن کاملPreservation of Adiabatic Invariants under
Symplectic methods, like the Verlet method, are a standard tool for the long term integration of Hamiltonian systems as they arise, for example , in molecular dynamics. One of the reasons for the popularity of symplectic methods is the conservation of energy over very long periods of time up to small uctuations that scale with the order of the method. In this paper, we discuss a qualitative fea...
متن کامل